1 



(N 
O 
O 
(N 

CD 
00 

(N 
(N 

(N 
> 

m 

I 

0^ 
O 
(N 
O 

(D 



X 



Temporal meson correlators at finite temperature on quenched 
anisotropic lattice* 
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"^Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

We study charmonium correlators at finite temperature in quenched anisotropic lattice QCD. The smearing 
technique is applied to enhance the low energy part of the correlator. We use two analysis procedures: the 
maximum entropy method for extraction of the spectral function without assuming specific form, as an estimate 
of the shape of spectral function, and the fit assuming typical forms as quantitative evaluation of the parameters 
associated to the forms. We find that at T ~ O.QTc the ground state peak has almost the same mass as at T^O and 
almost vanishing width. At T ~ l.lTc, our result suggests that the correlator still has nontrivial peak structure 
at almost the same position as below Tc with finite width. 



1. Introduction 

It is expected that the hadronic excitation 
modes strongly change their properties around 
the QCD phase transition. Such changes of 
hadron properties may signal the occurrence of 
phase transition in heavy ion collision experi- 
ments. With the potential model approach, it 
has been expected that the charmonium masses 
are shifted in the vicinity of Tc ■ In the plasma 
phase, they are expected to dissolve, and resul- 
tant suppression of J/tp formation has been re- 
garded as one of most important signals of plasma 
formation Q. Therefore it is important to verify 
these phenomena using lattice QCD in a model 
independent way. 

In principle, information of excitation modes 
can be extracted from the Matsubara Green func- 
tion in the Euclidean time direction. In practice, 
however, there are several problems at finite tem- 
perature. One problem is that the degrees of free- 
dom is restricted by the short temporal extent. 
To avoid this difficulty, we use an anisotropic lat- 
tice, on which temporal lattice spacing is finer 
than the spatial one. The other problem is that 
one need to extract the low energy structure from 
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the correlator at the region of 0(1/T), where the 
correlator contains contribution from wide range 
of frequency of the spectral function. We circum- 
vent this problem by applying the smearing tech- 
nique, which enhances the low energy part of the 
correlator. 

In order to extract a reliable information on the 
spectral function, we suggest to use maximum en- 
tropy method (MEM) D and method m a 
complementary manner. The former has an ad- 
vantage that it does not require to assume spe- 
cific form for the spectral function. Once the 
form of spectral function is estimated, however, 
the latter approach is more quantitative for eval- 
uation of properties of the mode, such as mass 
and width. We first discuss how the applicability 
of these procedures to finite temperatures is jus- 
tified using correlators at T = by varying num- 
bers of degrees of freedom used in analyses. Then 
we discuss changes of spectral function below and 
above Tc- Details of this work will be presented 
in future publication. 

2. Lattice setup 

We use quenched lattices of the sizes 20'^ x Nt, 
where Nt — 160, 32, and 26 which roughly cor- 
respond to T 2± 0, O.OTc, and l.lTc, respectively. 
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The zero temperature lattice, and the setup of 
lattice parameters are the same as those used in 
Ref. 0. The gauge configurations are generated 
with the Wilson plaquette action at /3 = 6.10 with 
the anisotropy = aa/a-t = 4. The spatial cutoff 
set by the hadronic radius ro is = 2.030(13) 
GeV. Nt = 28 is close to the phase transition. 
The quark action is the 0{a) improved Wilson ac- 
tion Q], with the hopping parameter correspond- 
ing to the quark mass itlq — 0.98 GeV. 
The correlator is represented as 



C{t) 



{O{x,t)O\0,0)). 



The operator 0{x,t) is 
Oix, t) = J2 't>iy)(i{x + y, t)Tq{x, t), 

V 



(1) 



(2) 



where 4x4 matrix T specifies the quantum num- 
ber and 0(y) is the smearing function for which 
we use the wave function observed at r = 0. The 
correlators are measured on 500 configurations at 
T = and 1000 configurations at T > 0. In order 
to reduce the statistical error, we average 16 cor- 
relators measured with different source points on 
each configuration. 

3. Analysis procedure 

We focus on the low energy structure of the 
spectral function. We use two types of analysis 
method: the maximum entropy method (MEM), 
and the standard fit with ansatz for the shape 
of spectral function A(u}). MEM has an advan- 
tage that it does not require to assume a specific 
form for the spectral function. When the num- 
ber of temporal points is sufficiently large, and 
the data region is also in sufficiently long dis- 
tance, MEM can reproduce the spectral function 
successfully However, when the number of 

temporal points is small, the extracted spectral 
function has a large uncertainty quantitatively, 
and sometimes even qualitatively. On the other 
hand, once a rough estimate of shape of the spec- 
tral function is in hand, fit gives more quanti- 
tative result for the parameters of assumed form. 
Therefore, complementary use of these methods 
are preferable for justification of their applicabil- 
ity and quantitative analysis. We use MEM for 
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Figure 1. The spectral function in the vector 
channel determined by MEM at r = 0. 

rough estimation of a shape of spectral function, 
and fit for more quantitative analysis. 

In MEM analysis we test several values of toq 
for the default model function, m{uj) = niQUi'^, 

= 1 is 



and fitting ranges [t 



where tr 



fixed. For the fit, we adopt two forms for each 
peak of spectral function: S function (denoted 
as pole) form, and the relativistic Brcit-Wigner 
{BW) form. 



AiLo) 



(3) 



Combining them, we apply 2-pole, 1-BW, and 
BW+pole fits to the correlators. In BW+pole fit 
the S function is used for subtraction of contri- 
bution from larger lo region (such as an excited 
state). 

4. Results of analysis 

In the following, we show only the result for the 
vector channel, while similar result is observed 
for the pseudoscalar channel. In order for MEM 
works at T > 0, it should produce at T = a 
structure which is stable under the change of fit- 
ting range, tmax- We regard this as a necessary 
condition for applicability of MEM to problems 
at T > 0. Figure shows the result of MEM 
at T = with several values of tmax- The po- 
sition of peak corresponding to the ground state 
is stable, while the width becomes broader for 
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Figure 2. Result at T > for the vector chan- 
nel. Top panel shows the result of 1-BW and 
BW+pole fits for the ground state peak width at 
T ~ 0.9Tc and l.lTc. Bottom panel is the result 
of MEM at T ~ l.lTc. 



width 7 consistent with zero. 

Now we turn to the study of spectral function 
at T > 0. At 0.9Tc, MEM gives a similar result 
as T = case. The BW+pole fit gives almost the 
same mass as at T = 0, and small value of width 
(top panel of Fig. Since this is the most gen- 
eral fit form, we conclude that at T ~ 0.9Tc, the 
ground state has almost the same mass as at T = 
and almost vanishing width. The other fits sup- 
port this conclusion. The 2-pole fit is rather well 
applied, while the obtained mass of the ground 
state slightly decreases as increasing tmin- The 
result of 1-BW fit approaches to the consistent 
value with BW+pole fit as tmin increases. 

The result of MEM at T ~ LIT^ is displayed 
in Fig. 2 (bottom). Although the width is larger 
than T < Tc, the result still exhibits a peak struc- 
ture around the same energy region as at T < Tc. 
Therefore we perform the same types of fit analy- 
sis as at T < Tc. The 2-pole fit and BW type fits 
give inconsistent results, and the latter fits indi- 
cate that the spectral function has a peak with 
almost the same mass at T < Tc and the width 
of order of 200 MeV, as shown in Fig. 2 (top). 
This result indicate that the charmonium corre- 
lator may still have a nontrivial structure slightly 
above Tc In this analysis, we use specific 

smearing function. In order to verify that the 
observed spectral function is not an artifact 0, 
comparative study with different smearing func- 
tion is necessary, which is now underway. 



smaller tmax- Although the position of the first 
excited state changes with tmax, we do not con- 
sider this seriously, since we are interested only in 
the ground state, and the contribution of excited 
states are small owing to the smearing of oper- 
ator. Therefore, using only restricted numbers 
of points, which are inevitable at T > 0, MEM 
can produce a consistent result with the case of 
sufficiently large tmax, at least for a rough esti- 
mate of lowest peak structure. This is in contrast 
to the case of a correlator of local operators, for 
which our MEM analysis does not produce a sta- 
ble result under the same change of tmax- The 
correlator is well fitted to the 2-pole form, and 
both 1-BW and BW+pole fits give the values of 
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